Astrophysical Journal, 536, June 10, 2000 



DYNAMICS OF WARM ABSORBING GAS 
IN SEYFERT GALAXIES: NGC 5548 

Mark C. Bottorff 

Department of Physics & Astronomy, University of Kentucky, Lexington, KY 40506-0055 

Kirk T. Korista 

Department of Physics, Western Michigan University, Kalamazoo, MI 49008 

and 

Isaac Shlosman 

Department of Physics & Astronomy, University of Kentucky, Lexington, KY 40506-0055 

ABSTRACT 

A hydromagnetic (MHD) wind from a clumpy molecular accretion disk is invoked to explain 
observations of warm absorbing gas at UV and X-ray energies in Seyfert 1 galaxies. This paper 
focuses on two important issues: (1) compatibility of kinematics and dynamics of an MHD wind 
with the observed properties of warm absorbers; and (2) the relationship between the UV and 
X-ray absorbing gases. We provide an in-depth comparison between the MHD wind model and 
the well-studied Seyfert 1 galaxy NGC 5548, which at high spectral resolution exhibits a number 
of discrete UV absorption components. Model parameters used in this study have been fixed 
by our previous work on the variability of broad-line region in this object. The detailed UVX 
absorption structure in NGC 5548 is modeled in order to infer the position, size, density and 
temperature of each component, as well as the overall dynamics of the absorbing gas. 

We find that for NGC 5548: (1) the total column densities of O vii, O viii and H, inferred 
from X-ray observations are reproduced by constraining the UV ion column densities of C IV 
and N V in each component to lie within a factor of 2 of their observed values and optimizing 
over the possible sets of component ionization states and C iv column densities; (2) the warm 
absorbing gas exists in the outer part of the wind and is not a continuation of the flow in the broad 
emission-line region; and (3) the warm absorber extends both in radial and polar directions and is 
ionization stratified. X-ray absorption is found to be heavily biased towards smaller radii and UV 
absorption originates at larger distances from the central continuum source. Our analysis shows 
that the discrete absorption components along the line-of-sight are intrinsically clumpy. Density 
differences between kinematic components result in a range of ionization and recombination time 
scales. 

We further test the applicability of the MHD wind to warm absorbers in general, by construct- 
ing a quasi-continuous flow model, and extending it to arbitrary aspect angles. Constraining the 
ionization and volume filling factor for a generic case, we estimate the fraction of Seyfert Is 
having detectable warm absorbers with larger O Vii column density than O Viii, and the range 
of total hydrogen column densities. We also find that the ratio of O Vli to O Viii optical depths 
can serve as a new diagnostic of AGN aspect angle. 

Finally, the thermal stability of the UVX absorption model is tested. We find that all kine- 
matic components in NGC 5548 are thermally stable to isobaric perturbations. In a more general 
case, we show that the magnetic field is crucial in order to stabilize the warm absorber gas over 
a wide range of incident continuum spectral energy distribution and gas metallicity. 
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1. INTRODUCTION 

Highly-ionized, tenuous, absorbing gas has been 
detected in the soft X-ray spectra of many active 
galactic nuclei (AGNs) at low and moderate {AE/E = 
0.02 — 0.1) spectral resolutions. At least half of the 
Seyfert 1 galaxies observed by the ASCA satellite 
indicate the presence of this "warm absorber" gas 
(hereafter WA) in their X-ray spectra (e.g., review by 
Mushotzky 1997). Other AGNs, such as radio-quiet 
and radio- loud quasars, host WA as well (e.g.. Green 
and Mathur 1996, Mathur and Elvis 1995, Sicbert et 
al. 1996, Ulrich et al. 1999), though no statistics are 
currently available. In addition, broad absorption- 
line quasars (hereafter BAL QSOs), which constitute 
at least 10% of all quasars, are conspicuously weak in 
the soft X-rays — a possible signature of a WA gas 
(Mathur, Elvis, & Singh 1995; Green & Mathur 1996; 
Gallagher et al. 1999). 

The X-ray WA gas is inferred from the softening 

of the continuTim at the low X-ray energies, between 
0.6 keV and a few keV, and is due to many unre- 
solved absorption edges from highly ionized species 
of elements, such as O, Ne, etc. The strongest edge 
depths correspond to O vii and O viii at 0.72 and 
0.85 keV, respectively. It is generally found that the 
total hydrogen column density of absorbing material 
lies between 10^^ cm^^ and lO^""* cm^^, and may even 
reach 10^^ cm^^ (Reynolds et al. 1997, hereafter R97; 
George et al. 1998, hereafter G98). The energy reso- 
lution of the ASCA satellite data does not allow for 
a study of the kinematics of the WA gas based on 
the X-ray absorption profiles. The thermodynamic 
state of the gas, however, including its temperature 
and density, can be inferred with a moderate amount 
of modeling and invoking observed X-ray variability 
(Otani et al. 1995, Reynolds & Fabian 1995; R97). 

Intrinsic blueshifted UV absorption was detected 
in AGNs by the lUE (Bromage et al. 1985; Koratkar 
et al. 1996; Leech et al. 1991; Wahcr et al. 1990; Voit 
et al. 1987; Veron et al. 1985) and HUT (e.g., Kriss 
et al. 1992). High and sometimes low- ionization lines 
are observed, e.g., C iv, N v, O vi. Si iv, Mg ii, 
and others. The higher S/N and spectral resolution 
data taken by the HST found many more cases of 



blueshifted intrinsic absorption in the UV resonance 
lines. Preliminary statistics indicate that intrinsic 
UV absorption is present in at least half of Seyfert 1 
galaxy spectra (Crenshaw et al. 1999 [C99]). In addi- 
tion to the BAL QSOs, there arc some quasars that 
exhibit intrinsic narrow absorption features (see, for 
example, Mathur, Wilkes, & Elvis 1999). 

In all known cases, the UV resonance line ab- 
sorption is present in sources with known soft X- 
ray absorption, hinting of a connection between the 
two types of absorption (C99; Mathur et al. 1999; 
Brandt, Laor, & Wills 1999), though the nature of 
this connection is not yet clear. The correspondence 
between the UV and X-ray absorbing gas was ad- 
vanced by Mathur et al. (1995), on the basis of single 
slab modeling of UVX absorbing gas in NGC 5548, 
using HST/FOS. ROSAT, and ASCA data in their 
analysis. Similar analysis was performed for other ob- 
jects (e.g., NGC 3516, Mathur et al. 1997; NGC 3783, 
Shields & Hamann 1997). The mean FOS spectrum 
of NGC 5548 represented the best UV data at the 
time, and its low spectral resolution (~ 200 km s"^) 
justified the single slab approximation. Much higher 
spectral resolution ffSTGHRS and STIS/Echelle have 
since revealed the presence of a number of distinct 
kinematic components within the outfiowing mate- 
rial, with FWHM of 50 kms'^ - 160 kms'^ (Cren- 
shaw & Kracmcr 1999). These high resolution spec- 
tra of NGC 5548 have uncovered at least five separate 
C IV absorption components (C99; Crenshaw & Krae- 
mer 1999; Mathur, Elvis, & Wilkes 1999 [MEW99]). 

The presence of a structure within the intrin- 
sic absorption adds another degree of complexity to 
the kinematical models of AGNs. It is important, 
therefore, to understand the relationship between the 
outflowing highly-ionized gas, the broad emission-line 
(BEL) and absorption-line regions in these objects. 
Numerous models have been proposed to address this 
problem (e.g., Blandford 1990; Arav, Shlosman & 
Weymann 1997; Konigl & Kartje 1994). In this paper 
we extend the hydromagnctic wind model of Emmer- 
ing, Blandford and Shlosman (1992, hereafter EBS), 
which provides a dynamical explanation for the the 
properties of the broad emission-line regions (BLRs) 
in AGNs (Bottorff et al. 1997, hereafter Paper 1), to 



-4- 



the UVX absorption phenomena. In the EBS model, 
gas in a dusty molecular accretion disk is loaded on 
magnetic field lines and centrifugally launched form- 
ing a magnetized wind. The basic difference between 
the EBS model and other outflow models is that it 
provides for a 3D wind dynamics, with helical mo- 
tion around the axis being superposed onto the radial 
motion. The wind is clumpy and stratified both in 
density and ionization. Paper 1 addressed the ra- 
dial wind structure at smaller latitudes above the 
optically-thick disk. First, it was found that the 
BLR is geometrically thick and has substantial opti- 
cal depth. As a result, the inner part of the wind be- 
comes progressively ionized with increasing latitude, 
while the outer part of the wind remains cold and 
molecular as long as it is shielded by the BLR. As 
the outer material continues to rise it becomes ion- 
ized and eventually passes the observer's line-of-sight, 
where the wind is detected via absorption in both UV 
and X-ray energy bands. The linc-of-sight velocity 
of the wind depends on its injection radius on the 
disk, so a range of velocities is expected. The hydro- 
magnetic wind model, therefore, naturally explains 
both the UV emission and absorption systems seen 
in C IV and other lines, as well as X-ray absorption. 
While our approach is general, we emphasize the re- 
sults within the context of observations of NGC 5548. 

In Section 2, we summarize the important obser- 
vations of the WA gas, including the column densi- 
ties of O VII, O VIII and total H from the X-ray data, 
and the column densities of C iv, N v and H i from 
the UV data. In Section 3, we consider specifically 
the case of NGC 5548 and its discrete UV absorption 
components. The velocity field of the EBS model 
adopted in Paper 1 is then used to determine the lo- 
cation and line-of-sight width of the UV absorption 
components. The observed column densities of N V 
and C IV in conjunction with the location and line-of- 
sight width of each component allow us to determine 
the density, temperature, volume filling fraction and 
column densities of O Vii, O VIII, H I and total H 
via photoionization modeling. The calculated ionic 
column densities are compared with observations, to 
determine whether or not the X-ray absorption in 
NGC 5548 is compatible with the model and in or- 



der to determine the association between the UV and 
X-ray warm absorbers. 

Section 4 describes a generic WA model. The 
multicomponent absorption, as seen in the C IV UV 
absorber of NGC 5548, is smoothed into an equiva- 
lent single absorbing column, to study the effects of 
a distributed WA in Seyfert galaxies without specific 
knowledge of kinematic absorption components. The 
generic model is investigated for arbitrary aspect an- 
gles to the observers. Wc estimate the fraction of 
AGNs having detectable WA gas, and larger O viii 
than O VII optical depth. Finally, in the Appendix, 
we discuss the general thermal stability of magnetized 
UVX absorbing gas. 

2. OBSERVATIONAL SUMMARY 

The spectral resolution of ASCA is insufficient 
to clearly resolve and measure oxygen and other ele- 
ment absorption edges without additional modeling. 
This is achieved by producing simplified photoioniza- 
tion models of WA gas using observed AGN continua. 
Free parameters in the model are subsequently tuned 
to obtain the best fit for the calculated transmitted 
spectrum to the observed one. 

A common approach taken by R97 and G98 makes 
a number of essential assumptions, namely, the WA 
is considered (1) to be a thin slab; (2) to have a con- 
stant total hydrogen number density; (3) to be in 
thermal equilibrium; (4) to be exposed to a power- 
law soft X-ray continuum of photon index F ~ 2, and 
(5) to possess solar abundances. The total hydrogen 
column density (or equivalently slab thickness) Nh, 
the X-ray ionization parameter Ux and T are then 
varied until an optimal fit to the data is obtained. 
Additional free parameters include Galactic absorp- 
tion, a fraction of unattenuated central continuum, a 
covering factor of the emitting gas, etc. 

The main inferred properties of WA gas from X- 
ray observations are the optical depths of O vii and 
O VIII, and the column density of total H along the 
line-of-sight. Ionization stratification in WA gas due 
to differences in density and location is indicated in 
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at least one object, namely MCG 6-30-15. The latter 
exhibits variability in the O VIII edge depth which 
is anticorrelated with the continuum, and little or no 
variability in O VII (Fabian et al. 1994, 1995; Otani et 
al. 1995). The difference in response is frequently at- 
tributed in literature to the recombination time scale, 
t{Xi), of ion Xi given by 

t{Xi) = l/a{Xi,T,)ne, (1) 

where rie is the electron density, and a{Xi,Te) is the 
recombination coefficient, which is a function of the 
ion species Xi and a slowly varying function of elec- 
tron temperature Te (but see section 3.3.5). In such 
a case, over a wide range in temperature a{0 vii, Tg) 
and a{0 Viii,Te) are of the same order of magni- 
tude, so recombination time differences between the 
two ions are mainly due to differences in density. The 
response of O Vlll and the lack of response of O vii in 
MCG 6-30-15 during the observing time is interpreted 
as evidence that O vii absorption occurs in a lower 
density gas than the O Viil absorption. Photoion- 
ization modeling and absorption depth variability in 
MCG 6-30-15 also indicate that the O vii and O viii 
absorption occurs at different locations (Otani et al. 
1995). An additional indication of ionization stratifi- 
cation in WA AGN gas is that objects with significant 
optical reddening exhibit largest O vii absorption — 
a hint that X-ray absorption may occur in regions 
where dust can survive (Komossa & Fink 1997). 

If the UV and X-ray absorbers are related, then 
some of the assumptions used in previous modeling 
efforts must be modified in view of the detection of 
multiple kinematic components in Seyfert galaxies, 
in particular in NGC 5548, which exhibits 5 sepa- 
rate, blueshifted UV absorption components. Thus, 
we abandon the single zone model in favor of a mul- 
tiple zone constrained by the EBS model. The con- 
stant density approximation is relaxed by estimating 
density for each of the kinematic components sepa- 
rately, by placing them at different distances from 
the central continuum source, based on kinematical 
constraints of the EBS model. We assume thermal 
and ionization equilibrium in the UVX gas (but see 
Krolik & Kriss 1995, and Nicastro et al. 1999). Tar- 
get values for the total hydrogen column density, Nh, 



and the optical depths t{0 vii) and t(0 viii), which 
we convert to column densities, are adopted from the 
X-ray observations discussed in R97 and are listed in 
Table la. 

In addition to the observations at X-ray energies, 
we also consider the HST UV GHRS observations of 
the 5 discrete, blueshifted absorption components of 
C IV, N V and H i detected in NGC 5548 on 1996 
Aug. 24 and on 1996 Feb. 17 (C99) (see Table lb). 
Wc adopt the ionic column densities measured from 
the GHRS observations over the STIS observations 
(Crenshaw & Kraemer 1999) because the former are 
closer in time to the ASCA X-ray observations of 
NGC 5548 on 1993 July 27 (R97). The single excep- 
tion is for the column density of H i which is lacking 
in component 1 in the GHRS spectra so we adopt the 
value from the STIS data. Note that the UV obser- 
vations are separated by seven months and the X-ray 
observations are separated from the UV observations 
by about 2 years. Extensive multiwavelength mon- 
itoring campaigns of this object (Clavel et al. 1991; 
Korista et al. 1995) show variations of UV emission 
lines which correlate with the continuum over times 
of days to weeks. C iv emission in NGC 5548 has, 
for example, varied in strength by a factor of ^ 2 in 
a month's time and C IV absorption equivalent width 
may be similarly variable (Shull & Sachs 1993). Cau- 
tion must therefore be taken when using noncontem- 
poraneous data to determine the ionization state of 
the absorbing gas. We also point out that reported 
values of the ionic column density depend upon the 
procedure used to fit absorption line profiles. For 
example, given the same set of HST/ GHRS spectra 
of the C IV absorption systems, MEW99 and C99 
derived values of the C IV column densities of com- 
ponents 1 — 3 which differed by 0.2 — 0.5 dex. Next, 
modeling that goes into determining observationally- 
derived ionic column densities does not explicitly ac- 
count for the possible effects of continuum scatter- 
ing into the absorption troughs, although these ef- 
fects do appear grossly in the estimate of the incom- 
plete line-of-sight continuum coverage. This model- 
ing must also make assumptions concerning the un- 
derlying emission from the narrow and BEL regions. 
Finally, the absorption lines are often observed to be 
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saturated, yet non-black. The line-of-sight covering 
fraction of the ions along the various line-of-sight ve- 
locities must be estimated from the degree of satu- 
ration present. This is usually done by comparing 
the depths of the resonance line doublet pairs. Dif- 
ferences in line-of-sight effective continuum coverage 
amongst the various ions can be expected and are 
often found. The H l column densities reported by 
C99 and Crenshaw & Kraemer (1999) suffer from a 
lack of constraints in this respect (Lya is a singlet, 
and no higher order Lyman line data are available), 
and should therefore be considered as lower bounds 
(Crenshaw, private communication). Thus the un- 
certainties in the ionic column densities reported in 
the literature are internal uncertainties only, the total 
uncertainties must be significantly larger. 

Our adopted ionic column density set is listed 
in Tables la and lb. Given the above considerations, 
we have disregarded their reported uncertainties, and 
instead assigned a factor of 2 uncertainty above and 
below their reported values. An exception is com- 
ponent 1 of C IV, which at the limit of detection is 
treated as a lower bound (Crenshaw, private commu- 
nication). These assigned ranges of uncertainty about 
the reported ionic column density values, which we 
will refer to as observational constraints, will be used 
to find a plausible set of solutions to the UVX ab- 
sorber in NGC 5548. 



3. DYNAMICS AND EMISSION 
MODELING FOR AGNs 

3.1. Hydromagnetic Wind Model 

The magnetohydrodynamic (MHD) solution for 
a stationary, axisymmetric flow in cylindrical coordi- 
nates is given by Blandford & Payne (1982) and de- 
scribes a self-similar, cold, nonrclativistic MHD flow 
from an idealized Keplerian disk around a point mass. 
This solution can be written in terms of variables %, 
^, ^ and a scaling parameter ro which are related to 
cylindrical coordinates via r = [?'o^(x), (j), rax]- Here, 
X is the coordinate along a field line, £,{x) is found as 
part of a self-consistent solution to the MHD equa- 



tions, and To is the field line footpoint on the disk. 
The flow velocity components are given by 

V = W{x)f{x),9{x),f{x)]^/GM/^o, (2) 

where a prime denotes differentiation with respect to 
X, and M is the mass of the central black hole. The 
EES provide for a non-trivial extension of Blandford 
& Payne solution for an arbitrary scaling of volume 
density n oc r^^ and magnetic field B oc Tq^'''^^^^^. 
At the base of the flow, the rotational velocity, v^, is 

— 1/2 

Keplerian so t;,^ oc rg . The functions C(x)i fix) 
and g{x) are chosen to satisfy the flow MHD equa- 
tions subject to the above scalings of p, B and v^. 
This implies that the Alfvcn speed scales with the 
disk Keplerian velocity, and the specific angular mo- 
mentum and energy in the flow will scale similarly 
to their Keplerian counterparts, while the disk mass 
loss per decade in radius scales as cx Tq ^'^\ The 
above velocity law is of course a non-Keplerian one 
at all finite distances from the disk (e.g., EES and 
Paper 1). 

EES utilized an analytic approximation for ^(x) 

in which ^(x) ^ X^^'^^ ^^nd asymptotically tends to 
the Blandford & Payne solution. To attach the solu- 
tion to the disk at ro, ^(x) was constructed to have 
the form ^(x) = Vx/c2 + 1, so that ^(0) = 1. The 
constant C2 = (tan 0o)/2, where is the initial angle 
between the magnetic line and the disk. The MHD 
equations have been solved for /(x) and ^(x)) so as 
to be consistent with the analytic form of ^(x)- The 
model parameters have been fixed in Paper 1 by mod- 
eling the BLR response to the variation of the central 
continuum (see section 3.2). We use b = 1.5, so n, the 
particle density, scales as n oc 

where R is the spherical radius from the central mass. 
We show below that this choice in b is consistent with 
the variation in density inferred from observations, al- 
though the number of observational points is small. 
This result, therefore, should be taken with necessary 
caution. We have further neglected the eff'ect of ra- 
diation pressure on the gas dynamics (e.g., de Kool 
& Begelman 1995), as it will not qialitatively change 
our results. 
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3.2. Constraining the MHD Flow with the 
BLR Spectrum 



With an analytical solution for the volume cmis- 
sivity EBS showed that realistic BLR emission line 
profiles can be produced. Paper 1 continued these ef- 
forts replacing the analytic emissivity function with 
an emissivity function obtained from a fit to an amal- 
gam of optically thick clouds calculated by the pho- 
toionization code Cloudy (Ferland 1996). Emission 
line anisotropy and finite optical depth effects were 
introduced to model C IV emission from optically 
thick BLR clouds in the well studied object, NGC 5548. below. 
As part of the modeling, variations in C IV emission in 
response to observed variations of the continuum, in- 
cluding light travel time effects, were studied. A time 
series of synthetic C IV emission line profiles were gen- 
erated and compared with the observed C iv emission 
line profiles from the HST (Korista et al. 1995), and a 
model fit to the data was thereby achieved. The BLR 
parameters in NGC 5548 deduced from the model 
were the central black hole mass (3 x IO'^Mq), physi- 
cal extent of the C iv emitting gas (1 to 24 light days) 
and the orientation of the observer line-of-sight rela- 
tive to the axis of symmetry (40°). Here, we extend 
the results of Paper 1 and analyze the flow beyond 
the BLR and at larger latitudes above the disk, in 
particular as the flow crosses the observer's line-of- 
sight. 



flow, which is the continuation of the BLR flow, as it 
crosses the linc-of-sight to the observer. Due to rather 
unfavorable conditions for loading and acceleration at 
the base of the wind, it is natural to expect a number 
of separate kinematic components along the linc-of- 
sight. The five blueshifted UV absorbers reported in 
C99 and MEW99 for NGC 5548, while perhaps not 
fully distinct from one another, are interpreted here 
as kinematically distinct regions containing gas in an 
EBS-type magnetized wind. Moreover, each individ- 
ual component is expected itself to be clumpy and is 
characterized by a volume filling factor, as we show 



A point on a streamline is defined in terms of 
the footpoint radius ro, the spherical radius R, and 
the observer's aspect angle i with respect to the disk 
axis, namely 

ro = Q{i,eo)R, (3) 



where 



Q{i,Oo) = 



\J cos^j ■ 



tan (?osin i — cosi 



tan^^n 



(4) 



In Paper 1, for the best fit model of the BLR in 
NGC 5548, the launch angle and aspect angle i, 
were fomid to be 20° and 40°, respectively. In this 
case, R « 10.4ro. Figure 1 illustrates the model ge- 
ometry for NGC 5548. A patch of gas filaments lifts 
off the disk between cylindrical radii r\ and r-i and 



rises along a helical trajectory, forming a parabola in 
In the ideal EBS flow, cold molecular gas is launched, , „ i .n 

_ _ . _ _ _ the rar-plane. Gas that starts at r\ eventually crosses 

the observer line-of-sight at a spherical radius R\ and 
gas that starts at r-i does this at R^. 



from a Keplerian disk and is flung along the magnetic 
fleld lines, like beads along a rotating wire. Most of 
the magnetic lines, however, are inclined unfavorably 
and hence will not be able to accelerate the gas. The 
centrifugally accelerated gas becomes illuminated by 
the central ionizing continuum at some latitude. The 
distribution of densities and ionization parameters, 
combined with the MHD model's velocity fleld pro- 
duces typically observed broad emission-line profiles. 
For NGC 5548, Paper 1 found that the BLR Hes 
within a toroidal wedge ±30° of the equatorial disk 
plane. The boundary of the wedge was interpreted 
as coming from optically-thick wind filaments experi- 
encing thick-to-thin transition. In this section, we fo- 
cus on the optically thin (to the ionizing continuum) 



The number density of a total hydrogen along 
the fiow line that starts at ro is prescribed by the 
model to be 



n 



riA 
m 



-1.5 



(5) 



Here ua is a normalization constant, setting the den- 
sity n on the innermost fiow line at m = 1.0, where m 
is the Alfven Mach number. The value of m increases 
from unity at about 17° above the disk to about 
lO'^'^ at 85° above the disk. The ionization parame- 
ter U (X 1/nR^, therefore, scales as t/ oc inQ'^{i,0o) 
along a flow line. At small latitudes, the effective ion- 
ization parameter for lower column density filaments 
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Fig. 1. — Schematic representation of the warm ab- 
sorber column. Gas, destined to contribute to the 
warm absorber, rises off a disk in the z = plane 
between n and r2 and crosses the observer's line-of- 
sight (OLS) between Ri and i?2, where it is detected 
via continuum absorption. Gas inside Ri is overion- 
ized, so there is little contribution to the absorption 
column densities of O VII, O VIII, C IV or N V in 
this region. Outside i?2 the flow has been cut off by 
mechanisms described in the text. 

will be even smaller, due to finite optical depth effects 
in the BLR (Paper f). The change in radius and den- 
sity causes the ionization parameter along any flow 
line to have a maximum at about 30° above the disk, 
where the BLR flow becomes optically thin. Here, 
the BLR flow is no longer self-shielded, becomes ove- 
rionized and remains so at larger latitudes, where it 
crosses the observer's line-of-sight. As a result, this 
gas does not contribute significantly to the UVX ab- 
sorbing column. The gas filaments launched at larger 
radii than the BLR fiow are less ionized and form the 
bulk of the UVX absorption. This absorbing gas is 
indicated by the thick line segment in Figure 1. 

While the inner boundary of the wind is set by 



photoionization, the outer boundary can be fixed by 
considering the amplification of magnetic fields in the 
disk. An efficient generation of magnetic fields in the 
disk will cause its expulsion through the buoyancy ef- 
fects and probable field reconnection above the disk 
(e.g., Galeev, Rosner & Vaiana 1979). The energy 
release during the reconnection can in principle in- 
ject disk material and load it onto the large-scale 
field treading the disk, resulting in the hydromagnctic 
wind. While the large-scale disk-threading poloidal 
field is probably governed by global considerations, 
the disk-generated fields are expected to be ampli- 
fied by differential rotation. Such dynamo effect will 
operate only within the "radius of infiuence," rsH, 
of the central black hole (BH), because for r > tbh 
the inner galactic rotation is typically that of a solid 
body. Hence, the outer wind cutoff can plausibly re- 
sult from a dramatic decrease in the shear outside 
Tbh, which prevents the magnetic field in the disk 
from reconnecting and loading the molecular mate- 
rial onto large-scale field lines. This radius can be 
estimated from equating the gravitational potentials 
of the BH and that of the host galaxy, namely, rsH ~ 
1.3 X 10"^^m7t;^^ cm, where m-j = A-Ibh/^O'^ and 
V2 = vk/WO kms~^ is the rotational velocity in the 
galactic disk. Adopting the BH mass in NGC 5548 
from Paper 1, leads to i?2 ~ 1-4 x 10^" cm for a typ- 
ical vr = 2.5. As we show below in § 3.3.1 this is 
within a factor of ~ 2 — 3 from the outer boundary 
of the observed kinematic components. 

3.3. Multi-Component Absorption Model 
for NGC 5548 

3.3.1. Location of Absorption Components 

The positions of separate absorption components 
are calculated using the relationship between R and 
Vobs, the line-of-sight component of wind velocity. 
The velocity Vobs is obtained by projecting the veloc- 
ity field onto the observer line-of-sight and is given 
by 

Vobs = K'(x)/(x)sim + /(x)cosi] VGMA^. (6) 
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Substituting ro = Q(40°, 20°)i? and Af = 3 x 10^ Mq 
for NGC 5548 (Paper 1), gives 



Assuming solar abundances, un = 9.33 x 10 ^ and 
ac — 3.55 X 10^^, we obtain the ion fraction ratio 



i? 7.3 X 10^^ 



/ Vobs \ ' 

U03 kms-iy 



(7) 



Observed component velocities in the rest frame of 
NGC 5548 from C99 and the FWHM of each compo- 
nent which we use below are listed in Table 3b. In- 
serting the observed component velocities into Equa- 
tion 7, yield component distances which range from 
about 2 to about 87 parsecs (Table 3b) . The most im- 
portant consequence is that the UVX gas in NGC 5548 
appears to be positioned well outside its BLR, both 
in radii and in latitude. Hence, within the framework 
of the EBS model, which neglects the effects of radi- 
ation pressure, the the UVX absorbing gas originates 
in the outer part of the wind, beyond the BLR flow. 



3.3.2. Calculation of Possible Ionization States 

The observational constraints on the ionic col- 
umn densities and their ratios in section 2 are now 
used to find the limits on the ratios of the N v ion 
fraction, f{N v), with the C iv ion fraction, /(C iv). 
These estimates allow the determination of a set of 
possible gas ionization states, i.e., the particle den- 
sity, electron temperature and various ion fractions 
in each of the UVX absorption components. We as- 
sume a total hydrogen number density n as given by 
equation 5, so the column density N{Xi) of ion Xi is 

NiX,) = / tf{X,)axndR « tf{X,)ax I ndR, 
Jar Jar 

(8) 

where e is the volume filling factor in the gas column 
which covers the central continuum source, f{Xi) is 
the ion fraction, ax is the abundance of element X, 
and Ai? is the line-of-sight geometrical width of the 
cloud. If the ionization parameter is a slowly varying 
function of radius with little change over Ai? and e is 
constant within each component, then ef{Xi)ax can 
be taken out of the integral. Then f{N v)/f{C iv) 
for each component is given by 



N{N v) 
N{C IV) 



f{N V) GN 

f{C IV) ac ' 



(9) 



fjN V) 
fiC IV) 



3.805 



N{N v) 
N{C IV) ■ 



(10) 



The central values of log[/(A^ v)/f{C iv)] for 
components 1 — 5, assuming the column densities 
listed in C99 and solar abundances, are listed in the 
final column of Table 1. In accord with section 2, 
we assign factors of two uncertainty to their values. 
Component 1 is a special case because the reported 
values for N{C iv) provide a lower bound only. 

The range of possible ionization states for each 
component was determined using Cloudy (version 90.04). 
The luminosity and spectral energy distribution (SED) 
of the continuum source in NGC 5548, combined with 
the distance of each component (Eq. 7), is used to 
calculate the incident flux ^{H) of hydrogen ionizing 
photons falling on each absorption component. We 
used the following continuum shape for the SED of 
NGC 5548, i.e.. 



and 



U{2 keV) 



403.3° 



(11) 



(12) 



/^(2500 A) 

where kTjfi = 0.01 Ryd is the spectrum cutoff in 
the near-infrared, and a in Equation 11 is chosen 
to satisfy Equation 12.|^ The input parameters are 
log {Tbb) = 6.683, aox = -1.20, auv = -1-20 and 
ax = -0.90, chosen to be in accord with lUE/ROSAT 
observations of the NGC 5548 continuum (Walter et 
al. 1994). Using the above SED, the ionizing luminos- 
ity of the continuum source is estimated at lO*'* '^ erg s~^ 
{Hq = 75 kms~^ Mpc~^), based upon the mean con- 
tinuum flux at 1350A, obtained during the 1993 ob- 
serving campaign (Korista et al. 1995). 

With the above assumptions, the ionization pa- 
rameter U (= ^[H]/nc), which controls the ioniza- 
tion state of the gas, is tuned by adjusting the hy- 
drogen number density n in a series of Cloudy cal- 
culations for an optically-thin slab of width 10^°cm. 



1 Outside the range 1.36 eV-100 keV, a = 0. Above 100 keV, we 
assume fu oc . 
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in this manner, a set of synthetic photoionization gas 
states covering a range of possible /(TV v)//(C iv) 
values is produced for each component. 

In order to provide for a complete description of 
the UVX absorption in each kinematic component, 
it is necessary to know the column densities of dif- 
ferent ions. Note that, if wc replace N{N v) with 
N{Xi) and un with ax in Equation 9, we obtain 
a generic expression for the column density of ion 
Xi. In the approximation of Equation 9, the quan- 
tity [f{Xi)/f{C iy)]{ax/ac) is constant, so N{Xi) oc 
N{C IV). Thus for fixed f {N v)//(C iv), all ionic 
column densities scale directly with the C iv col- 
umn density. The parameters f{N v)/f{C iv) and 
N{C iv) then map out the set of possible ionic col- 
umn densities for each component. We use these two 
parameters to determine the electron temperature Tg 
and column densities of N v, O vi, O vii, O viii, H i 
and H, on a grid of [f(N v)/f(C iv),iV(C iv)] pairs 
that are consistent with the observational constraints. 

The grid of solutions constrained by f(N v) /f(C iv) 
and N(C iv) were used to determine the range of pos- 
sible values of the column densities of O vi, O vii, 
and O VIII in each component. The factor of 2 un- 
certainty in the observational constraints on the C IV 
and N v column densities and their ratios translates 
into column density uncertainties of a factor of ~ 10 
on O VI, a factor of ~ 100 on O vii and a factor of 
~ 1,000 on O Vlll. The inferred ranges on the col- 
umn densities of H and H i are no better. H i ranges 
over a factor of ~ 100 and H over a factor of ~ 1, 000 
(even more in component 1). We, therefore, obtain an 
optimal solution for UVX absorption in NGC 5548, 
within the five component sets of possible solutions. 
Since individual column densities of C IV and N V for 
each component are observed, but only total column 
densities are available for O VII and O VIII, we search 
for a solution that minimizes the chi-squares, Xo' ^'^^ 
oxygen ion column densities, namely 

2 (E -=1^(0 VII),- -iV(OVIl)o6.)2 

Xo = — -2 

"O VII 

, (E.tiMO vni)i-jV(Qvni)o6«)^ 

"O VIII 



Here N{0 vii)j and N{0 Viii)j are the calculated col- 
umn densities of O Vil and O Vlll in a component j 
(see Table 3), N{0 vii)obs and N{0 viii)obs are the 
observed total column densities of O vii and O viii 
(see Table 1), and ctq vn and ao vm are the adopted 
observational uncertainties (section 2) in the O Vii 
and O viii column densities. Our search of the solu- 
tion grid yields optimal Xo — 0.13. The total column 
density of O vii is overpredicted by 30% and the to- 
tal column density of O Viii is underpredicted by 9%, 
compared to target values in Table la, placing both 
well within the observational constraints. 



3.3.3. An Optimal Solution for UVX Absorption 
using the EBS Model 

The model column densities of C iv, N v, O vi, 
O vii, O viii, H I and H in each component and their 

sum over all components are listed in Table 3a. Com- 
ponent C IV and N v column densities automatically 
satisfy the observational constraints because the op- 
timization was done on a grid of values bounded by 
the constraints. The component column density val- 
ues for O VI, O VII, O VIII, and H in each component 
are predictions of our model. The good match be- 
tween the predicted and observed total O vii, O viil, 
and hydrogen column densities, given present obser- 
vational constraints, indicates that the soft X-ray WA 
and the UV resonance-line absorber are directly re- 
lated. The EBS model constrains the bulk of the WA 
to lie within two of the kinematic "components" of 
the outflow, separated by about 500 km s~^ in line- 
of-sight velocity. While observations at the required 
resolution do not yet exist to test these predictions in 
detail, the recently launched X-ray satellite Chandra 
may be able to detect the O vii and O viii absorption 
edges of the individual kinematic absorption compo- 
nents in NGC 5548, and the recently launched FUSE 
satellite should at least detect the presence of absorp- 
tion in O VI. 

The predicted column densities of H i from our 

optimal solution do not match well their reported 
values (Table lb). The observed component column 
density values are factors of 4 — 28 below their coun- 
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terparts in the optimal solution. However, as men- 
tioned in section 2, the obscrvationally inferred H I 
column densities are the least certain, and are lower 
limits to their true values, at best. 

The hydrogen number density for each compo- 
nent in the optimal solution as well as Tg, U and other 
physical values are listed in Table 3b. The density de- 
rived for each component using Cloudy depends on 
the kinematic structure of the EBS model via Equa- 
tion 7, but is independent of the parameter 6 = 1.5 
used to describe the scaling of the EBS density in 
Equation 5. It is therefore possible to test whether 
the assumed value of b used in Paper 1 and here in 
this paper is appropriate. 



^^^^ 


J log(n) 




log(U) 





IB. 5 19 19.5 20 20.5 



log(R[om]) 

Fig. 2. — Logarithmic plot of density and ionization 
parameter as a function of _R[cm] for the discrete com- 
ponent model of NGC 5548. The straight lines are the 
least squares fits. 

Equations 3-5 together imply that along any ray 
through the origin n oc R~^'^, and, therefore, the 
ionization parameter varies as U on R^'^'^. If indi- 
vidual absorption components obey this density pre- 
scription, then n and U from component to compo- 
nent should vary accordingly. As a test, we fit log(n) 
and \og(U) from the optimal solution at the center 
of each absorption component, as functions of log(i?) 
with a linear least squares fit. Figure 2 shows log(n), 
log(L'') and the least squares fit of both. The range in 
uncertainty is ±0.2 and corresponds to our photoion- 



ization grid's resolution in gas density. We find that 
the linear least squares fits to log(n) and log{U) vs. 
log(i?) have slopes of -1.57 ± 0.16 and -0.45 ± 0.16, 
respectively. These are consistent with our specific 
EBS assignment of n oc R~^-^ and U oc R~^-^. Since 
(7 is a slowly decreasing function of R, our simpli- 
fication of the integral in Equation 8 is reasonable, 
assuming that there are no significant optical depth 
effects along the absorbing column and that in each 
component AR/R is sufficiently small, so that there 
are no significant changes in U over the width of a 
component either. The picture that emerges for the 
wind model of NGC 5548 is that of optieally thin fil- 
aments originating over a large range of radii in the 
disk, from 0.2 pc to about 9 pc. These filaments be- 
ing centrifugally accelerated along the magnetic field 
lines, cross the line-of-sight about ten times further 
out. Along the line-of-sight, the ionization parameter 
drops as a weak function of distance [U oc 
but because the components arc spread out over a 
factor of ~ 40 in radius, there is a factor of ~ 6 drop 
in ionization parameter so the wind is ionization- 
stratified along the line-of-sight. The thermal stabil- 
ity of kinematic components in NGC 5548 is analyzed 
in the Appendix. 

3.3.4. Absorption Component Widths and Volume 
Filling Factors 

A more complete physical picture of the UVX ab- 
sorption system is provided by determining the line- 
of-sight width, Aii, of each absorption component. 
In this section we estimate the ratio AR/R, discuss 
the radial gas distribution, separation of individual 
components, their degree of dumpiness, variation of 
U, etc. 

The line-of-sight extension of each kinematic com- 
ponent can be estimated from the FWHM of the ab- 
sorption components (given in C99) and the electron 
temperature, provided by the optimal solution (Ta- 
ble 3). The FWHM is attributed to the effects of 
thermal Dopplcr broadening and velocity gradients 
in the MHD flow. The component due to velocity 
gradients, Av, is obtained by making a correction to 
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the FWHM by subtracting the thermal broadening 
which we take to be Vg = ^2kTe/mi, and by taking 
the ion mass nii to be the average mass of carbon 
and nitrogen atoms. This average is taken because 
the FWHMs Hstcd in C99 were the result of averag- 
ing the FWHM determined from the N V and C IV 
absorption lines. Table 3b lists Vg and Av, defined by 
At; = FWHM—Vs- All of the Av values are positive, 
therefore, some of the line width cannot be attributed 
solely to thermal Doppler broadening. Equation 7 is 
used to determine the spatial extent AR of each com- 
ponent via 



AR = R 



AR 



Vobs 



Av 
~2~ 
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Vobs 
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Fig. 3. — The location and width of the discrete ab- 
sorption components of the optimal model for warm 
absorbing gas in NGC 5548. 

Numerical values for AR in each component arc 
also listed in Table 3b. Figure 3 shows the width 
and location of each component for the optimal so- 
lution plotted in log(ii). Gaps between the compo- 
nents indicate that they are well separated (except 
for possibly components 3 and 4), and that AR val- 
ues grow monotonically with R. For the first three 
components AR/R < 0.32, but it is of order unity for 



components 4 and 5. According to the EBS model, 
the ionization parameter will vary across the compo- 
nent by Ufront/Uback = ^ {R + ARl2) / {R - AR/2). 
The values of R and AR for each component (Ta- 
ble 3b) yield percentage differences in U measured at 
the front and back sides of each component. These 
are ~ 11%, 6%, 17%, 71% and 61% for components 1 
through 5, respectively. Thus our approximation in 
Equations 8, 9 and 10 is reasonable only for the first 
three components. For components 4 and 5 we esti- 
mate that the effect will result in an error of about 
—0.1 dex in the N v and C IV column densities, and 
as much as -1-0.4 dex uncertainty in the oxygen col- 
umn densities. Fortunately, these errors are minor 
in those components where oxygen ions predominate 
(i.e., 1 and 3), and, therefore, the breakdown of our 
assumption of constant U across individual compo- 
nents does not affect much the predicted integrated 
oxygen column densities. 

An estimate for the volume filling factor is ob- 
tained from the column density. Using n oc R~^'^ 
as a density profile, the integral in Equation 8 is 
given by J^f^ndR = n{R)2R[{l- AR/2R)-'^^'^ - {1 + 
AR/2R)~^/'^], where n{R) and R are the density and 
radius at the center of each component of width AR. 
Hence 

N{Xi) « 

tf{Xi)axn{R)2R[{l-AR/2R)-^/^-{l+AR/2R)-^''^] 

(15) 

If n{Xi) is the density of ion Xi at i?, then 



n{Xi) = f{Xi)axn{R). 



(16) 



Substituting Equation 16 into Equation 15, letting 
q = AR/2R and solving for e gives 



N{Xi) 



9(1 



j2)l/2 



n{Xi)AR (1 + g)V2 _ (1 _ qy/2 ■ 



(17) 



The quantity N{Xi) /n{Xi)AR is the filling fac- 
tor, when the density is constant throughout the com- 
ponent. The term involving q makes a correction for 
the variation in density across a component. The 
correction factors, however, are close to unity. For 
components 1 through 5, they are 0.994, 0.998, 0.984, 
0.844 and 0.871, respectively. This requires the choice 



- 13 - 



of an ion to define e, and we choose C IV since it is 
already used as one of the solution parameters. We 
note that N v or even a blend of the C IV and N v 
column densities could be used as well, but the dif- 
ferences in calculated e are not significant. 

Calculated values of log(e) for the optimal so- 
lution range from —5.58 to —1.37 (Table 3b). This 
means that an absorption component is clumpy and 
has a filamentary structure. In other words, each of 
the five distinct absorption components in NGC 5548 
consists of filaments which move along similar trajec- 
tories and have probably a common origin in the disk. 
On the other hand, we presume that the distribution 
of points of origin of any given absorption component 
in the disk plane is random. The dumpiness of indi- 
vidual components combined with the geometry and 
non-radial nature of the MHD flow has an important 
implication for the longevity of UV absorption strTic- 
tures in AGNs, which we discuss in the next section. 

3.3.5. Temporal Changes in the UVX Absorbing 

Column 

Additional information about the character of 
the absorbing column of gas is inferred from the tem- 
poral behavior (or lack thereof) of UV and X-ray 
absorption. Changes in absorption lines in princi- 
ple can be due to the motion of material into or out 
of the observer's line-of-sight, or it may be due to 
changes in the ionization structure of the gas brought 
about by changes in the incident continuum or both. 
The X-ray absorption line variability has been dis- 
cussed in Reynolds & Fabian (1995), Reynolds (1997) 
and Otani et al. (1995), where ASCA observations 
of MCG 6-30-15 show apparent anti-correlation of 
the O VIII edge depth with an increase in the con- 
tinuum level. A similar conclusion was reached by 
Reynolds & Fabian (1995). Observations of UV ab- 
sorption lines in AGNs also reveal that they change 
with time. Shuh & Sachs (1993) used lUE data to 
show that the C IV absorption equivalent width in 
NGC 5548 responds to continuum changes in an anti- 
correlated fashion on a time scale less than about 4 
days (this value is likely to be rather uncertain due 



to the fact that the temporal sampling of this mon- 
itoring campaign was four days). UV observations 
of other AGNs spaced over many months reveal that 
absorption lines change in depth with time, but the 
relationship with the continuum is not clear. While 
the depth of an absorption line may vary, there are no 
unambiguous observations of a change in the velocity 
centroid of the line (Weymann 1997), which argues 
against purely radial motion, since one would expect 
clouds moving only radially to accelerate along the 
line-of-sight, resulting in time with a shift in the ve- 
locity centroid. This is in agreement with our mod- 
eling of gas dynamics in the BLR of NGC 5548 (Pa- 
per 1) using the EBS wind, which ruled out a purely 
radial motion as well. Here we focus on the temporal 
phenomena associated with UV and X-ray absorption 
in AGNs. 

Our analysis of the discrete absorption compo- 
nents in the previous section suggests that each com- 
ponent is an agglomeration of filaments moving along 
the same trajectory. The heuristic model is that the 
wind filaments are confined by the ambient magnetic 
field and are extended along the field lines. The ex- 
act details of gas loading onto the magnetic lines are 
subject to future work. Filaments which rotate out 
of the line-of-sight are replace by their neighbors, so 
the velocity centroid of the absorption component is 
not expected to change with time. The time scale 
tc over which we expect an absorption structure to 
survive is given by tc = AR/v^, where AR is the 
width and is the rotational velocity of the absorp- 
tion component. For the first UV absorption compo- 
nent in NGC 5548, where tc has the smallest value, 
we find ~ lO'* km s^^ and AR ~ 10^^ cm, lead- 
ing to tc ^ 3, 000 years. The proposed model can, 
therefore, explain the observed long constancy of the 
velocity ccntroids of UV absorption structures. 

Given the long dynamical time scale, observed 
temporal changes in the absorbing column are at- 
tributed instead to changes in the photoionization 
structure of the absorbing gas, rather than a change 
in gas density along the line-of-sight, and are presum- 
ably due to fluctuations in the incident continuum. 
The response of the absorption within the EBS wind 
to continuum fluctuations is investigated by consider- 
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ing the effect of a sudden drop in the continuum level 
on the column density of C IV, O vii and O VIII. C IV 
is chosen to compare with observation (Shull & Sachs 
1993), and O vii and O viii are chosen to compare 
and contrast the response of our model of NGC 5548 
with observations of MCG 6-30-15. 

In the above scenario, ions will begin to recom- 
bine with the drop in the continuum intensity. How- 
ever, the time scale for recombination as given by 
Equation 1 docs not account simTiltanc;ously for the 
cascade into the population of Xi ions from the pop- 
ulation of Xi+i ions, and the cascade out of the pop- 
ulation of Xi ions into the population of ions. 
Equation 1 is hence misleading and is replaced for 
consistency by 



t{Xi) = 



1 



a{Xi)n( 



fix.) 



a(Xi. 



(18) 



The value of t{Xi) in Eq. 18 establishes the m,ini- 
mum time, that gas will respond to a decrease in the 
continuum, and depends on the local density and on 
the ion population. As a single measure within a 
component, we take the average ion density-weighted 
recombination time which is 



< t{Xi) >= 



_ J^j^t{Xi)nXiedR 



N- 



Xi 



(19) 



This definition weights more heavily recombination 
times in regions with higher column density contri- 
butions. Assuming the validity of Equations 8—10, 
nXi,nXi+i and rie are all proportional to i?^^-^. In 
addition, within each component is approximately 
independent of R and, therefore, a{Xi) and a{Xi-i) 
are also independent of B. The result is that t{Xi)nx, 
is independent of R and may be factored out of the 
integral. This gives 



< t{Xi) >= 



tjXj) \RnxdR)eAR 
Nxi 



(20) 



where t{Xi) \ji and njf. (i?) arc t{Xi) and njf. , re- 
spectively, evaluated at R. Hence, < t{Xi) > may be 
positive or negative, depending on whether or not the 
column density of ion X^ is anti-correlated or corre- 
lated with a decrease in the ionizing continuum, re- 
spectively. 



Calculated values of < t{Xi) > for C iv, O vil 
and O VIII in NGC 5548 are listed in Table 4. Not sur- 
prisingly, these recombination time scales differ from 
those calculated from Equation 1 — an issue some- 
times ignored in the literature. For a fixed electron 
density, the ionic distributions as well as the effect of 
recombinations to lower ionization states can produce 
a wide range of the recombination time scales for the 
various metallic ions. In Table 4, the computed value 
of < t{Xi) > for each ion differs from component to 
component. The recombination time scale for C IV 
spans a factor of ~ 5, 000. The ranges for O vii and 
O VIII, span factors of ^ 600 and ^ 50, respectively. 
Whether or not changes of ion Xi in a single compo- 
nent have a significant effect on the total ionic column 
density over all components of ion Xi , depends on the 
relative contribution of Xi in the single component. 
(All reported variations in ionic UVX column densi- 
ties were measured over the total value, due to insuffi- 
cient spectral resolution.) For this reason the percent 
contributions of C IV, O Vli and O Viii are also shown 
in Table 4. The values are used to define the column 
density-weighted mean recombination time over all 
components, << t{Xi) >> (Table 4). The individ- 
ual recombination times in each component serve as 
lower bound response time scales to changes in the 
ionizing continuum flux, and << t{C iv) >> may be 
compared in a crude manner with observations. In 
our model, << t{C iv) >>« +3.7 days is reasonably 
consistent with the results of Shull & Sachs (1993), 
namely, that C iv absorption equivalent width is anti- 
correlated with continuum variations with a delay of 
roughly four days. 

Long term X-ray observations of NGC 5548 do 
not yet exist, so the time scales reported in Table 4 
for O VII and O viii cannot be tested. Instead wc 
compare the results of this model to the observations 
of MCG 6-30-15. In our model of the UVX absorp- 
tion flow in NGC 5548, the bulk of the contribution 
to the column density of these two oxygen ions re- 
sides in components 1 and 3. However, the values of 
<< t{0 vii) >> and << t{0 Vlll) >> are consid- 
erably skewed by the extremely long recombination 
time scale of component 5. As a conservative mea- 
sure, we, therefore, consider only components 1 and 
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3. In these components, the column density contri- 
butions are roughly equal for both O vii and O VIII, 
but the recombination times in the first component 
are over 10 times shorter. These recombination times, 
which represent the minimum possible response time 
scale for a continuum drop, are < t{0 vil) >« —3 
days and < t{0 viii) -50 days. (The O vii and 
O VIII edge optical depths are correlated with the 
ionizing continuum intensity.) We contrast this with 
observations of MCG 6-30-15 (see section 2). During 
the second half of the ~ 4 day observation period of 
this object, the X-ray continuum dropped by a factor 
of 2. The response was that the O viii optical depth 
anti-correlated with the continuum on a time scale 
of 10^ s, but there was no observed change in the 
O VII optical depth. 

A possible explanation for the observed behavior 
of the O VIII edge optical depth in MCG 6-30-15 is 
that the absorbing gas in MCG 6-30-15 is more highly 
ionized than in NGC 5548. The observed continuum 
in MCG 6-30-15 is much harder than in NGC 5548, 
and the WA gas may, therefore, support a large pop- 
ulation of O IX ions. Equation 18 allows then for a 
complex dependency of the recombination time scales, 
which might explain the differences between the pre- 
dicted behavor of NGC 5548 and the observed be- 
havior of MCG 6-30-15. Long looks by Chandra and 
other advanced X-ray telescopes should establish the 
relationships between the continuum and WA optical 
depths in AGNs. 



4. A GENERIC MODEL FOR WARM 
ABSORBING GAS 

The discrete UVX absorber model for NGC 5548 
does not allow for its direct application to other ob- 
jects. This is mainly because detailed kinematic data 
is not available, the model does not provide for the 
UVX absorption dependence on the aspect angle, and 
because it does not specify the appropriate boundary 
conditions not along the line-of-sight. To generalize 
the model, we proceed by spatially smoothing the 
contributions from individual UVX kinematic com- 
ponents of our model and by calculating the optical 



properties of the MHD flow along a number of aspect 
angles. The resulting generic model is used to esti- 
mate orientation effects on the soft X-ray absorption 
properties of AGNs. The general aspects of thermal 
stability of the WA gas in this model are discussed in 
the Appendix. 

4.1. Column Densities and physical 
conditions 




log(R cm) log(R cm) 

Fig. 4. — Left: Plot of the ion fraction as a function 
of radius for the continuously distributed model of 
NGC 5548. For clarity, C iv and N iv values have 
been multiplied by a factor of 10. Hi has been mul- 
tiplied by a factor of 10''. Right: Logarithmic plot 
of the particle density (cm~^) and temperature (°K) 
as a function of radius (cm) for the continuously dis- 
tributed model of NGC 5548. 

We first build a quasi-continuous model of MHD 
flow in NGC 5548 by smoothing the kinematic com- 
ponents along the line-of-sight, by imposing obser- 
vational constraints on the total column density of 
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hydrogen, and by retaining the (x R~^/'^ density pro- 
file from the discrete component model. For sim- 
plicity, the outer edge of the continuous flow was set 
at logi? w 20.6 which coincides with the outer edge 
of the fifth component. Two free parameters, i.e., 
the volume filling factor and the inner radius, which 
we used as Cloudy input parameters, were tuned to 
fit observations of the column densities of N{C iv), 
N{N v), Nh, N{0 vii) and A^(Ovin). A reason- 
able fit to the ionic column densities is found for 
log(iZi) w 17.7 and e ~ 10"^ (Table 2). The first plot 
in Figure 4 clearly illustrates ionization stratification, 
with the ion fraction O viii achieving a maximum at 
smaller radii, O Vii at intermediate radii, and both 
C IV and N v peaking near the outer edge. The sec- 
ond plot in Figure 4 shows the density and electron 
temperature as a function of distance. Note the tem- 
peratures range from ~ 1.6 x lO^K at larger radii to 
2 X lO^K at smaller radii, where the bulk of the 
WA resides. 



4.2. Orientation Effects in AGNs and Soft 
X-Ray Absorption 

We now investigate the orientation effects in the 

continuous wind model of NGC 5548 by varying the 
viewer aspect angle. Figure 5 shows the column den- 
sities of Af(^ v), N(C IV), N{H), N{0 vii) and A^(0 
as a function of viewer aspect angle i. The den- 
sity profile along each aspect angle was adjusted so 
that the MHD flow is normalized to the solution for 
NGC 5548 (Table 2), which corresponds to the as- 
pect angle i = 40°. Absorbing column densities seen 
by an observer with aspect angle i > 60° are within 
the angular wedge shadowed by the BLR, i.e., by hy- 
drogen column densities in excess of 10^^ cm~^. Al- 
though we restrict our discussion to angles at 60° 
or below, and to smaller column densities compati- 
ble with those in NGC 5548, larger column densities 
cannot be excluded for other objects. We adopt the 
generic model for NGC 5548 as a template, when ap- 
plied to other Seyfert 1 galaxies. 

Figure 5 shows that N(0 vii) and N{0 viii) de- 
crease with decreasing i, making detection of the oxy- 
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Fig. 5. Plot of the log of the column densities 
(cm^^) of N V, C IV, H, O vii and O viil as a fimc- 
MSB^i of observer aspect angle i for the continuously 
distributed model. Values with angles greater than 
60° do not include contributions to the column from 
optically thick BLR clouds. The H column is there- 
fore most likely to be an order of magnitude (or more) 
higher when i > 60°. 

gen edges more difficult for small i. Using the instru- 
ment error bar in R97, and assuming a 2a detection 
criteria, absorption edges of O vii and O viii will 
be lost for aspect angles less than i 30° with the 
disk axis. Since no objects will be seen for i > 60°, 
due to obscuration, we restrict our analysis to angles 
i < 60°. In this case, if all Seyferts are taken as ran- 
domly oriented templates of NGC 5548, then about 
~ 73% of them will have detectable WA gas. 

Figure 5 also shows that N{0 vii) drops faster 
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than N{0 viii) when viewed at smaller aspect an- 
gles, and drops below N{0 viii) for angles less than 
50°. The fraction of Seyfert 1 galaxies with detectable 
WA showing N{0 viii) > N{0 vii) can be estimated 
from the ratio of the solid angle subtended on a sphere 
between the detection limit 30° and 50°, and the solid 
angle between 30° and 60°. This ratio is about 61%, 
for current instrument sensitivity. For comparison, 8 
out of 16 Seyfert 1 galaxies listed in R97 (50%) have 
N{0 viii) > N{0 vii). (We used only Seyfert Is 
from R97 that have definite error bounds on both the 
tabulated values of optical depth tq vii and tq vm-) 
Finally, we note that the predicted range of Nh varies 
from ~ 10^^ cm~^ to ~ 10'^'^ cm~^ within the ob- 
servable range of aspect angles, and this compares 
favorably with the inferred range of Nh in R97. Of 
course, one can question this approach, when one ob- 
ject is used as a template to describe the whole pop- 
ulation of Seyferts, and the Seyfert sample in R97 is 
small. The model's consistency with the observations 
in these respects is nevertheless interesting. Taken at 
face value, the ratio of optical depths of O vii and 
O viii ions can serve as a new diagnostic for AGN 
aspect angles. 

4.3. Diffuse Emission from the UVX 
Absorber 

The UVX absorbing gas reprocesses and scatters 
the incident continuum into diffuse emission across 
the UV and soft X-rays. We used our "generic" UVX 
absorber for NGC 5548 to estimate the magnitude 
of this omission. We computed Cloudy models for 
thermal local line widths and one, in which micro- 
turbulence broadened the local line width [oturb = 
100 km s~^). In the latter scenario, optically thick 
lines are desaturated, and the cross-section for pho- 
ton pumping of the line's upper level (continuum res- 
onance line scattering) is elevated, enhancing the line 
intensities. Photon pumping may be especially im- 
portant for the soft X-ray lines because of their small 
Boltzmann factors in photoionized gas. However, even 
under the assumption that the observer can see the 
gas represented by our generic UVX absorber dis- 
tributed over 47r steradians over the continuum source. 



the diffuse emission is relatively small. This is ex- 
pected, given the reported small optical depths of 
the absorber in this object. We find that this gas 
contributes at most 20% to the observed UV nar- 
row hnc fluxes of Lya, C iv, and N v in NGC 5548 
(Goad & Koratkar 1998). Emission of this signifi- 
cance may be sufficient to partially fill in some of 
the absorption troughs. Depending upon the cover- 
ing fraction and the contribution from photon pump- 
ing the diffuse soft X-ray line emission could be also 
be significant in comparison to the absorption edge 
depths (Netzer 1993). We will have to await the re- 
sults from the upcoming missions of Chandra, XMM, 
and Constellation- X to get hard observational con- 
straints from this portion of the spectrum. 

5. SUMMARY 

We have applied the EBS hydromagnetic wind 
model from a clumpy molecular accretion disk to the 
problem of UVX-ray absorption systems observed in 
Seyfert galaxies. We have studied the state of the 
UVX absorbing gas, with a particular emphasis on 
NGC 5548, where five discrete UV absorption sys- 
tems seen in this object were modeled. Extending 
our previous work on the dynamics of the BLR gas 
in NGC 5548 (Paper 1), we first inferred the location 
of each component relative to the continuum source. 
By estimating the ratio of the ion fractions of N V 
to C IV from the observed column densities and the 
use of photoionization modeling, we have determined 
possible values for the density, electron temperature 
and the ion fractions of C iv, N v, O vi, O vii O viii, 
H and H i. Optimization over the set of possible val- 
ues yielded the ion column densities in each kinematic 
component and allowed us to test the radial and po- 
lar density distributions of the WA in NGC 5548, 
within the EBS framework. We have also estimated 
the size and volume filling factor of each component, 
using the FWHMs and electron temperatures, and 
found that the flow in each component is clumpy. X- 
ray absorption by O Vii and O VIII column densities 
seen in NGC 5548 have been accounted for in the five 
UV absorption components, though main contribu- 
tion to the X-ray absorption comes from two compo- 
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nents. We, therefore, were able to explain the UVX 
absorption columns within the framework of the same 
dynamical model. An additional important point is 
that the model parameters used here are the best 
fitting parameters from Paper 1 to explain the broad 
emission-line variability in NGC 5548 during the 1989 
and 1993 observing campaigns. 

Secondly, we find that the WA gas in NGC 5548 
lies at larger radii from the central source than the 
BLR gas and at larger altitudes above the disk. This 
means that WA exists in the outer parts of the hydro- 
magnetic disk wind, and as such is not a continuation 
of the BLR flow. As a result we do not expect it to 
contribute to absorption along the line of sight unless 
the aspect angle is near the obscuring torus. With re- 
gard to the spatial extent of the WA gas we find it 
extends in both radial and polar directions from the 
central continuum source, and is ionization stratified. 

Thirdly, we have modeled the WA in NGC 5548 
also as a continuous flow, in order to investigate the 
generic properties of the UVX absorption outflows 
in Seyfert galaxies. We find that orientation effects 
may fully account for the fraction of objects with de- 
tectable WAs. Furthermore, the model predicts that 
the ratio of optical depths of O viii to O vii can serve 
as a new diagnostic of orientation of an AGN, namely 
of its rotation axis — an issue of a particular interest 
in AGN theory. In addition, we find that the UV line 
emission from the UVX flow probably accounts for 
no more than ^ 20% of the NLR emission in lines, 
such as Lya, C IV, N v, and O vi, though diffuse 
emission in the soft X-rays may be signiflcant under 
some conditions. 

A thermal stability analysis has been carried out 
for the model. We find that all of the five components 
in NGC 5548 are stable, though two components (1 
and 3) lie near an unstable, to isobaric perturbations, 
region on the S-curve. By considering the effects of 
magnetic field on the thermal stability of the WA, we 
find that these components can be further thermally 
stabilized by a modest field, even in the presence of 
substantial continuum fiuctuations, even when con- 
sidering extremes in SED or metallicity. Such stabi- 
lization by magnetic fields may explain the ubiquity 



of WA gas in AGNs. 

As more objects are observed with more sensi- 
tive instruments (e.g., Chandra and XMM coupled 
with the HST), better statistical tests of the wind 
model can be carried out, involving determination of 
the fraction of AGNs with detectable WA gas, the fre- 
quency of objects in which the O vii column density is 
larger than the O viii column density, and a determi- 
nation of the range of Nh- Future higher resolution 
X-ray spectra (ASCA, XMM) will provide additional 
kinematic information, such as oxygen edge and res- 
onance line absorption velocities and help to verify 
the link between X-ray and UV absorption. 
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A. THERMAL STABILITY OF WARM 
ABSORBING GAS 

The WA gas is found in a wide range of den- 
sities, temperatures, ionization and metallicity. It 
seems plausible, that this ubiquity of WAs in AGNs 

is related to their thermal stability. Permeating mag- 
netic fields which are expected to play an important 
dynamical role in AGNs, can affect the thermal sta- 
bility of the absorbing gas. Therefore, along with 
Field (1965) and EBS, we discuss the general aspects 
of thermal stability of magnetized gas and its appli- 
cation to NGC 5548 and to a generic model of the 
WA, subject to AGN radiation field. 

A.l. The "S-curve" and Thermal Stability: 
General Considerations 

For a gas in thermal equilibrium, the heating rate 
per unit volume, G{n,T), is balanced by the cooling 
rate per unit volume, A(n, T). If is the net heating 
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Fig. 6. — S-curves for various parameters, (a) S- 
curves for the SED of NGC 5548 for densities 10^ 
10^, 10^ and 10^ cm ^ . The thick hne contains both 
the 10^ cm^-^ and 10^ cm^'^ S-curves. Note that 
above 10^ K all the S-curves merge, illustrating that 
the S-curve is relatively insensitive to differences in 
density, (b) S-curves corresponding to different inci- 
dent continuum SEDs. The S-curve labeled 1 is for 
the SED of MCG 6-30-15; the S-curve labeled 2 is for 
the SED of NGC 5548; the S-curve labeled 3 is for 
the Mathews & Ferland (1987) composite spectrum 
of quasars. The hardness of the SED decreases from 
left to right, (c) S-curves for the SED of NGC 5548 
with three different metal abundances. The solid line 
corresponds to solar abundances; the dashed line is 
for a metal abundance of 0.5 x solar; the dotted line is 
for an abundance of 2. Ox solar, (d) Location of the 5 
discrete absorption components observed by C99 on 
the S-curve for NGC 5548. Segments are offset for 
clarity. Components 1 and 3 lie close to an unstable 
(to isobaric perturbations) portion of the S-curve. 



H= G{n,T) -A{n,T), 



(Al) 



then equilibrium occurs when H = 0. If gas is slightly 
perturbed, i.e., H ^ 0, the gas will either heat up or 
cool down. The change in the gas temperature de- 
fines whether it is thermally stable (Field 1965; Kro- 
lik, McKee & Tarter 1981). The equation H ^ 
is more conveniently treated when mapped to the 
log T — log U/T plane where the equilibrium curve 
has roughly an S-shape. Thermal stability of gas de- 
pends on the details of the S-curve, and the path that 
a perturbation is allowed to move along. Figure 6a 
displays a series of S-curves corresponding to differ- 
ent hydrogen gas densities for the SED of NGC 5548. 
The curves are relatively insensitive to the gas den- 
sities considered and become virtually indistinguish- 
able above T « 10^ K. 

To analyze thermal stability of gas, one must be 
able to follow the path of the perturbed gas relative 
to the S-curve. Figure 7 shows a small segment of the 
S-curve for NGC 5548, having a negative slope in the 
vicinity of point P. To the right of the S-curve, the 
net heating H is positive, i.e., perturbed gas always 
heats up. To the left of the S-curve, the net heating H 
is negative, meaning that perturbed gas always cools 
down. Crossing the S-curve at the equilibrium point 
P are three paths (G, L and I), representing differ- 
ent alternatives for an evolving perturbation off the 
point P. Path I, a vertical line, represents an isobaric 
perturbation]^ G and L represent perturbations con- 
strained to paths with slopes greater (less negative) 
and with slope less (more negative) than the slope 
of the S-curve at P, respectively. Along I and L, a 
small perturbation off the S-curve leads to a dramatic 
change in temperature. Perturbations constrained to 
move along G, on the other hand, are thermally stable 
since they are driven back to P. Hence, the thermal 
stability of the gas depends on the slope of develop- 
ing perturbation in log T — log U/T plane. The gas is 
thermally stable when this slope is greater (less nega- 
tive) than the slope of the S-curve at P. For an ionized 



If luminosity and radius arc fixed tlrcn U/T oc 1/nT, tlicrcforc, 
U/T — constant implies nT is constant. Hence, the gas moving on 
vertical lines in the T vs. U/T plane evolves isobarically. 
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Fig. 7. — Schematic depiction of an unstable (to iso- 
baric perturbations) portion of the S-curvc. To the 
left of the S-curve cooling exceeds heating. To the 
right of the S-curve heating exceeds cooling. Differ- 
ent perturbation constraints (paths G, L, and I) are 
shown crossing the point P. Arrows show the direction 
gas will evolve along a constraint when perturbed off 
the point P. Path I corresponds to a thermally unsta- 
ble isobaric perturbation at P. Path L has a slope less 
than the S-curve at point P, therefore, perturbations 
constrained on this path are unstable. Path G has a 
slope greater than the S-curve at P, so perturbations 
constrained to this path are stable. 

gas permeated by a magnetic field, as shown below, 

perturbations move along inclined paths in this di- 
agram, e.g., as G or L in Figure 7. Since stability 
depends on the relative steepness of the S-curve, as 
compared to that of the the perturbation path, we 
first investigate the factors which influence the shape 
and, therefore, the slope of the S-curve. 



A.2. Factors Affecting the S-curve: SED 
and Metallicity 

The overall shape of the S-curve is affected by the 
SED. Figure 6b, shows the S-curvcs for three different 
SEDs with labels corresponding to SED with photon 
index F = 1.8 and energy cutoffs between 13.6 eV and 
40 keV (this SED was utiHzed in Reynolds 1995), to 
the SED used in our modeling of NGC 5548, and to 
the typical quasar SED deduced by Matthews & Fer- 
land (1987). The sequence is in terms of decreasing 
hardness of the continuum and the greatest effect is 
on the upper branch of the S-curve. Here the temper- 
ature is mainly controlled by Compton heating and 
cooling. In the limit of large log(?7/T), the S-curve 
is asymptotic to the Compton temperature Tq which 
is given Tq = h < v > /Ak, where < v > is the 
flux-averaged photon frequency. Higher < y > cor- 
responds to harder SED, and, therefore, higher Tq, 
explaining the sequence of decreasing upper branch 
temperatures. 

The horizontal separation between S-curves (Fig. 6b) 
is due to decrease in the continuum hardness, when 
the percentage of high energy photons {hv > 1 Ryd) 
that heat the gas per unit frequency interval decreases. 
Decreasing the hardness, allows the gas to cool down, 
unless compensated for by increasing the number of 
high energy photons, which is the same as increasing 
the ionization parameter. Thus, for a given value of 
T, the S-curve lies further to the right (i.e., it has 
higher U/T values which correspond to higher values 
of U) as the hardness decreases. In spite of a wide 
range of hardness, the slopes of the S-curves vary 
little, particularly in regions where the slope is nega- 
tive. We further note that in all cases considered here, 
the S-curve slopes is always less than —1, in regions 
where the slope is negative. This is an important 
point because, as we show below, the path that gas 
perturbation treaded by a magnetic field follows, has 
slope close to —1, resulting in thermal stabilization. 

The bumps and wiggles in an S-curve are due to 
the metallicity of the gas. If AGNs are fueled with 
galactic ISM, its composition may reflect intense lo- 
calized star formation which is typically concentrated 
within nuclear and circumnuclear regions and, there- 
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fore, uniformity of composition is not guaranteed. In 
addition, when grains are directly exposed to an in- 
tense central continuum, their mantles will be de- 
stroyed and only the graphite cores will be able to 
survive photo-desorption (Draine & Salpeter 1979) 
at distances and in the radiation field relevant for 
WA in NGC 5548. This introduces the possibihty 
of metallicity gradients within the WA environment. 
Gas abundances are important, because at WA tem- 
peratures line radiation by heavier elements, such 
as oxygen, constitute the dominant coolant. Thus 
metallicity can influence the thermal structure and 
stability of the WA gas. This is demonstrated below. 

Figure 6c shows S-curves for the SED of NGC 5548 
for gas with three different metallicities. At temper- 
atures above T k, lO^X, higher metallicity curves 
are shifted further to the left than lower metallicity 
curves in the log T vs. log U /T plot. This is because 
the presence of more metal ions gives the gas a net 
higher energy absorption cross section than that at a 
lower metallicity, so that a particular temperature is 
achieved with a lower ionization parameter. In addi- 
tion, extra electrons provided by metal enhancement 
are able to recombine with iJ+ (and iJe++), thus 
the effective ionization of the gas is lower for a given 
ionization parameter. For example, along a line of 
log T « 5.85, the ionization parameter increases from 
roughly 0.85 to 1.14 to 1.25, as the metallicity drops. 
Below T K, 10^ K, the order is reversed because the 
metal ions recombine and become more efficient cool- 
ers via line emission than heaters via absorption. In- 
spection of the Figure 6c reveals that the slopes of the 
S-curves for different metallicities are similar, and in 
the regions where the slopes are negative, they all 
are less than —1. So, metallicity differences, even as 
large as a factor of 2 greater or less than solar, will 
not affect our conclusions about thermal stability of 
magnetically-dominated gas. 



A. 3. Thermal Stability of Magnetized Gas 

For isobaric perturbations, there are no dynam- 
ical effects and the resulting evolution of the pertur- 
bation is controlled by thermal effects only. In the 



absence of a magnetic field, an isobaric path corre- 
sponds to a vertical line (path I in Figure 7). For 
this kind of perturbation, regions where the slope of 
the S-curve is negative are thermally unstable. The 
location of the discrete UV/X-ray WA model is shown 
in the log T vs. log U/T plane (Fig. 6d). The tem- 
perature range of absorption components 1 and 3 lies 
close to a region which is thermally unstable to iso- 
baric perturbations, so changes in continuum lumi- 
nosity may drive these components into an unstable 
region.^ Components 2, 4 and 5 are in a region where 
the S-curve slope is positive, so these components are 
thermally stable. 

The mere existence of kinematic components 1 
and 3 in NGC 5548 hints that the gas is thermally 
stable, despite being close to unstable region. Such 
a stabilization mechanism can be provided by mag- 
netic fields permeating the gas (Field 1965; EBS). In 
the presence of magnetic field, isobaric perturbations 
are constrained to evolve along the curve defined by 
nkT + i?^/87r = A where A is a constant. To ana- 
lyze thermal stability, we rewrite this expression in 
terms of T and and show that when plotted in the 
log T vs. log U/T plane, the slope is approximately 
equal to —1. Assuming magnetic flux freezing, the 
magnetic field scales with the gas density = c/i^n^, 
where </> is a constant. Substitution into the expres- 
sion for total pressure gives nkT + (f)^n'^ /Sn — A. We 
eliminate n by solving 



ion 



U ^ L_ 

T < Eion > 47ri?2c7lT 

for n. This gives 

E 



(A2) 



(A3) 



where E = Lion/ < Eion > ^itB?c. Substitution into 
the expression for total pressure gives 



1^2 



Ek 



f 87rT2(|)2 



A, 



(A4) 



The temperature uncertainty for each component is estimated by 
searching for the largest and smallest values temperature values in the 
(n,N(C iv)) solution grid that are adjacent to the optimal solution. 
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and solving for T yields 



(t>E 
\/8¥ 



A{^Y-Ek{^) 



To thermally stabilize the gas requires 



d log(r) 

d log(C//T) 



>5, 



(A5) 



(A6) 



where S, the slope of the S-curve in the log T vs. 
log U/T plane, and the derivative on the left hand 
side of Equation A6 arc evaluated at the point where 
Equation A5 crosses the S-curve, yielding 



(flog T 
d\og (I) 



= -1- 



(A7) 



where (3 is the ratio of gas to magnetic field pres- 
sures. Note, that when the magnetic field B — > 0, 
then fi ^ CO, and dlog T/dlog {^:) — oo, recover- 
ing the strictly isobaric case. At the other extreme 
lies the limiting case of the cold EBS flow, in which 
the magnetic pressure completely dominates the gas 
pressure (/? « 1), so dlog T/cflog (^) w -1. 

The minimum value of B required to stabilize the 
gas is found by (equating the slope of the S-curve to 
Equation A7, and solving for Pmin- Figure 8 shows 
^/l^min values as a function of temperature for the 
two regions of instability, for the SED of NGC 5548. 
For stable regions, no magnetic field is needed, so 
l//3mm 0. Components 1 and 3 are close to the 
thermally unstable temperature corresponding to the 
left peak in Figure 8. The maximum l//3mm occurs 
at log T = 5.55, where l/Pmin = 0.395, giving 



1 



Ml 

-T7^ = 0-395, 

nkT 



(A8) 



where B and n are the values for log T = 5.55. Thus, 
according to Equation A8, if the value of the magnetic 
pressure is about 40% of the gas pressure or larger the 
gas is stabilized. Cloudy calculations and our wind 
model yield a density of log n = 4.99 at the center 
of component 1. With this density a lower bound on 
the magnetic field required for stabilization can be 
estimated. Solving Equation A8 for B gives 




B = VStt X 0.395nfcT. 



(A9) 



log(T) 



Fig. 8.— vs. T for the S-curve of NGC 5548. /3 
is the minimum value required to stabilize gas at the 
corresponding temperature. 

In the regime of the left peak, the gas is stable at 
all possible temperatures, ii B 7.0 x 10""^ Gauss. 
It should be noted, however, that the bound on B 
is determined only from the stability requirement. A 
more stringent requirement from the EBS model is 
that the flow is a "cold" MHD flow, which means 
1/(3 » 1, superseding the above criteria. The im- 
plication is that the slope given by Equation A7 is 
nearly equal to —1. We have noted above that varia- 
tions in SED or metallicity do not significantly change 
the slopes of the negative parts of the S-curve, and 
that the slopes are always less than —1, correspond- 
ing to case G of Figure 7, in which the slope of the 
perturbation is larger than the slope of the S-curve. 
Thus, in the case of NGC 5548, as described by the 
EBS flow, regions of the S-curve that correspond to 
thermally unstable regions to isobaric perturbations 
are thermally stabilized by the presence of a magnetic 
field. 
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Table la 

ASCA R97 X-ray Observational Summary NGC 5548 
(Adopted Target Values'') 



^Og{NH) 


21.71 


\og{No vn) 


17.95 


log (iVo VIII ) 


18.20 



Table lb 

HST C99 UV Observational Summary NGC 5548 
(Adopted Target Values") 



comp. \og[N{N v)] log[iV(C iv)] \og[N{H i)] log[' 



1 


14.30 < 


13.04 < 


14.13*^ 


2.14 > 


2 


13.78 


13.45 


13.88 


0.91 


3 


14.59 


13.68 


14.25 


1.49 


4 


14.81 


14.46 


14.52 


0.93 


5 


14.04 


13.61 


13.11 


1.01 




ITj.lT 


14. (i2 


14.87 





" All values have assigned error bars of ±log(2.0) except for those with a <, which indicates a lower 
bound or >, which indicates an upper bound. 

^STIS observation (Crenshaw & Kraemer 1999) 
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Table 2 

Model Solution for the Continuously Distributed 
Warm Absorber in NGC 5548 



parameter 


log(value)" 


■f^min 


17.75 


r^max 


20.59 




6.455 


I^max 


2.253 


N(H) 


21.76 


N(H i) 


15.94 


N(C IV) 


14.70 


N(N v) 


15.13 


N(0 vii) 


17.94 


N(0 viii) 


18.22 


e 


-2.716 



" values in cgs units. 
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Table 3a 

Optimal EBS Solution Column Densities' 



comp. 


C IV 


N V 


VI 


VII 


VIII 


H 


H I 


1 


13.50 


14.30 


16.15 


17.73 


17.84 


21.29 


15.57 


2 


13.45 


13.48 


14.57 


14.51 


13.00 


18.14 


14.11 


3 


13.62 


14.41 


16.25 


17.80 


17.88 


21.33 


15.66 


4 


14.46 


14.51 


15.65 


15.65 


14.20 


19.22 


15.14 


5 


13.49 


13.74 


15.11 


15.56 


14.55 


18.88 


14.35 




14.62 


14.94 


16.58 


18.07 


18.16 


21.61 


16.00 



'^Logarithims of column densities in units of cm 



Table 3b 

Optimal EBS Solution Physical Parameters'* 



comp. 


^obs 


FWHM'' 




Av 


f{N V) 

f(C IV) 


R 


AR 


n 


T 


U 


e 


1 


-1060 


114 


9.8 


104.2 


1.39 


18.8 


18.1 


5.0 


4.87 


-0.03 


-1.8 


2 


-655 


45 


5.4 


39.5 


0.61 


19.2 


18.3 


5.4 


4.37 


-1.29 


-5.6 


3 


-518 


90 


9.6 


80.4 


1.37 


19.4 


18.9 


3.8 


4.86 


-0.05 


-1.4 


4 


-344 


156 


5.4 


150.6 


0.63 


19.8 


19.8 


4.2 


4.36 


-1.23 


-4.8 


5 


-165 


74 


5.9 


68.7 


0.83 


20.4 


20.4 


2.6 


4.44 


-0.85 


-4.1 



"Velocities are given in km s ^. All other entries are given as logaritlim of tlie value expressed in cgs. 
''Values are from C99. 
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Table 4 

Component Ion Recombination Times 
and Percentage Total Column Contributions 
for NGC 5548 



n 


C IV" 




VII 


% 


VIII 


% 


1 


2.71 


7.59 


5.39 


45.71 


6.62 


47.86 


2 


4.23 


6.76 


5.46 


0.03 


5.54 


> 0.01 


3 


3.96 


10.00 


6.60 


53.70 


7.80 


52.48 


4 


5.30 


69.18 


6.59 


0.38 


6.57 


0.01 


5 


6.39 


7.41 


8.14 


0.31 


8.30 


0.02 


Y." 


5.51 




6.43 




7.54 





" Values are the logt(s) of the absolute value of the mean component recombination time. All C iv 
responses are anticorrelated and all O vii and O viii responses are correlated with a continuum drop. 

^ Percentage of the total ionic column density. 

Values are the logt(s) of the absolute value of the mean recombination time averaged over all compo- 
nents and weighted by their relative column density contributions. 



